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Abstract. We investigate the dependence of composition, spectrum and angular 
distributions of ultra-high energy cosmic rays above 10 19 eV from individual sources on 
their magnetization. We find that, especially for sources within a few megaparsecs from 
the observer, observable spectra and composition are severely modified if the source is 
surrounded by fields of ~ 10~ 7 G on scales of a few megaparsecs. Low energy particles 
diffuse over larger distances during their energy loss time. This leads to considerable 
hardening of the spectrum up to the energy where the loss distance becomes comparable 
to the source distance. Magnetized sources thus have very important consequences for 
observations, even if cosmic rays arrive within a few degrees from the source direction. 
At the same time, details in spectra and chemical composition may be intrinsically 
unpredictable because they depend on the unknown magnetic field structure. If 
primaries are predominantly nuclei of atomic mass A accelerated up to a maximum 
energy E max with spectra not much softer than E~ 2 , secondary protons from photo- 
disintegration can produce a conspicuous peak in the spectrum at energy ~ E max /A. 
A related feature appears in the average mass dependence on energy. 
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1. Introduction 

The origin of ultra-high energy cosmic rays (UHECR) above 10 19 eV (= lOEeV) is a 
mystery since many years PjJ |2] . Several next-generation experiments, most notably the 
Pierre Auger experiment now under construction jH] and the EUSO project [I] are now 
trying to solve this mystery. 

Although statistically meaningful information about the UHECR energy spectrum 
and arrival direction distribution has been accumulated, no conclusive picture for the 
nature and distribution of the sources emerges naturally from the data. There is on the 
one hand the approximate isotropic arrival direction distribution 5\ which indicates that 
we are observing a large number of weak or distant sources. On the other hand, there 
are also indications which point more towards a small number of local and therefore 
bright sources, especially at the highest energies: First, the AGASA ground array 
claims statistically significant multi-plets of events from the same directions within 
a few degrees [HI E|, although this is controversial [Zj and has not been seen so far 
by the fluorescence experiment HiRes The spectrum of this clustered component 
is oc E^ 1 ' 8 and thus much harder than the total spectrum [Hj. Second, nucleons 
above ~ 70 EeV suffer heavy energy losses due to photo-pion production on the cosmic 
microwave background — the Greisen-Zatsepin-Kuzmin (GZK) effect [9J — which limits 
the distance to possible sources to less than ~ 100 Mpc jTO] - Heavy nuclei at these 
energies are photo-disintegrated in the cosmic microwave background within a few 
Mpc [H]. For a uniform source distribution this would predict a "GZK cutoff", a 
drop in the spectrum. However, the existence of this "cutoff" is not established yet 
from the observations |12j . 

The picture is further complicated by the likely presence of large scale extra-galactic 
magnetic fields (EGMF) that will lead to deflection of any charged UHECR component. 
Magnetic fields are omnipresent in the Universe, but their true origin is still unclear jI3] . 
Magnetic fields in galaxies are observed with typical strengths of a few micro Gauss, 
but there are also some indications for fields correlated with larger structures such as 
galaxy clusters Magnetic fields as strong as ~ \\iG in sheets and filaments of 

the large scale galaxy distribution, such as in our Local Supercluster, are compatible 
with existing upper limits on Faraday rotation [TH ITK] . It is also possible that fossil 
cocoons of former radio galaxies, so called radio ghosts, contribute significantly to the 
isotropization of UHECR arrival directions ^B] • 

Only recently attempts have been made to simulate UHECR propagation in 
a realistically structured and magnetized universe based on large scale structure 
simulations [T7| IT%j . These two simulations used different models for the EGMF: 
Whereas seed fields were continuously injected at shocks in Ref. [T7], Ref. [TH] started 
from uniform seed fields. In both cases, the seed field strength had to be normalized 
to reproduce observed rotation measures. Interestingly, the EGMF in Ref. jTTj is 
considerably more extended and leads to much larger average deflections than in 
Ref. PH] • The main reason is probably due to the different EGMF models. This suggests 
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that the influence of EGMF on UHECR propagation is currently hard to quantify. The 
two works agree, however, at least in the following: First, the observer is most likely 
situated in a region with relatively weak EGMF, < 10~ 9 G. In Ref. [T7J which used an 
unconstrained large scale structure simulation, this was deduced from comparison of 
predicted and observed UHECR anisotropics, whereas in Ref. JB] this followed from 
their constrained simulation which is supposed to lead to an approximate model of our 
local extra-galactic environment. Second, the fields are not strong and extended enough 
to allow an interpretation of the large scale isotropy of UHECR arrival directions by only 
one nearby source, as suggested in Ref. [JJJ]. Finally, both these works were restricted 
to nucleons. Heavy nuclei may further complicate the issue and lead to large deflection 
even in the EGMF scenario of Ref. [TB] . 

In the present paper we study UHECR spectra, composition and angular 
distribution from individual discrete sources injecting protons or nuclei into a highly 
structured EGMF which is significant only in the neighborhood of the source. This 
is different from the more idealized Monte Carlo simulations for individual sources in 
Refs. |2H1 [2JJ 122 123 121 in the following points: The EGMF in these earlier simulations 
was idealized by a Kolmogorov spectrum with a possible coherent component, with a 
field energy density either following a Gaussian profile (2211211 or being constant [201123, 
or as organized in spatial cells with a given coherence length and a strength depending 
as a power law on the local density [21]. Furthermore, the observer was immersed in 
fields of similar strength as the source, and these simulations were restricted to nucleons. 

In the context of a highly structured EGMF it was already found that the auto- 
correlation function at degree scales decreases with source magnetization and can thus 
be used as a signature for magnetic fields [25] ITT] . In the present paper we will find 
that even if off-sets of arrival directions from the source direction are moderate, spectra 
and composition are considerably modified if the source is surrounded by EGMF of the 
order 10~ 7 G on scales of a few Mpc. 

In Sect. 2 we describe our simulations, in Sect. 3 the results are presented, and we 
conclude in Sect. 4. 

2. Simulations 

We will specifically study two cases where the location of the observer and source are 
chosen from the large scale structure simulation used in Ref. [T7j, see Ref. [21] • In the 
first case the source is at a distance of 3.3 Mpc, similar to the starburst galaxy M82 
or the radio galaxy Centaurus A. The field strength at the observer is ~ 3 x 10~ n G, 
whereas the source is immersed in fields of ~ 1CT 7 G. We will denote this scenario 
by "M82". In the second case, the source is immersed in a galaxy cluster at 18.6 
Mpc distance, reminiscent of the Virgo cluster. The field strength at the observer is 
~ 10~ n G, whereas the source is immersed in fields of a few 10~ 7 G on scales of several 
Mpc. We will denote this scenario by "Virgo". These two cases are visualized in Fig. ^ 
We consider the steady state situation which amounts to assuming that these sources are 
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active on time scales at least as long as the propagation time. In this case observables 
are time independent. 
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Figure 1. Log-scale two-dimensional cuts through magnetic field total strength (color 
scale, in Gauss) with source position indicated by the asterisk in the center. Left panel: 
Scenario "M82" . The observer is represented by the white disk to the left and above 
the source. Right panel: Scenario "Virgo" . The observer is represented by the grey 
disk in the upper right corner. 

The trajectory simulations are performed as described in Ref. |T2J. These 
simulations were generalized to include heavy nuclei and follow all secondary 
trajectories produced by photo-disintegration reactions, in the way described in 
Ref. [2H]. Interactions taken into account thus include pair production by protons, 
pion production, and photo-disintegration on the combined cosmic microwave, 
infrared/optical and radio photon backgrounds. For the latter, standard estimates were 
used, see Ref. [2H] for details. 

For simplicity, we will assume that nuclei are accelerated with a spectrum oc E~ 2 
up to a fixed maximal energy -E max = 4 x 10 21 eV, and we will restrict ourselves to 
either proton or iron primaries at injection. These are reasonable assumptions for radio 
galaxies [27] • We note that active galaxies are likely to accelerate heavy nuclei to ultra- 
high energies, see Refs [T7J EHj for discussions. 

The observer is modeled as a sphere and a trajectory is accepted and observationally 
relevant quantities are saved every time it crosses this sphere. For each scenario 10 6 
arriving trajectories were simulated. Galactic magnetic fields lead to deflections of 
maximally a few tens of degrees above 10 19 eV and can therefore be neglected as long 
as one is mainly interested in spectra and chemical composition. 

3. Results 

Fig. El compares the solid angle integrated spectra predicted by scenario "M82" for iron 
primaries, with and without EGMF. The left panel corresponds to the source position 
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Figure 2. Steady state all-particle spectra predicted by scenario "M82" for iron 
primaries with injection spectrum oc _E~ 2 - up to 4 x 10 21 eV. The blue curves are for the 
EGMF surrounding the source shown in Fig.^ left panel, and the red curve is without 
EGMF. Shown for comparison are the solid angle integrated AGASA [2U (dots) and 
HiResT j3U] (stars) data. The solid straight line marks the injection spectrum. The 
right panel is for a source off-set by about 1 Mpc from the position shown in Fig. ^ 
left panel. All fluxes have been normalized at 60 EeV. 
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Figure 3. Same as Fig. |2 left panel, but for proton primaries. Here the fluxes have 
been normalized at 300 EeV. 

shown in Fig. ^ left panel, whereas the right panel is for a slightly different source 
position for comparison. Fig. El shows the same for proton primaries for the source 
position shown in Fig. Q left panel. The normalizations are chosen such that this 
individual source can contribute at least partly to the measured UHECR flux at the 
highest energies. At lower energies cosmological sources are likely to dominate the 
observed flux [T7j . 

There are several notable features for the nearby source "M82" which we will discuss 
in turn: It is first of all apparent that whereas without fields the spectrum is only 
slightly modified due to the small source distance, the magnetized case leads to a strong 
modification of the spectrum. This can be understood by comparing the time delay due 
to propagation in the EGMF, shown in Fig. EJ with typical energy loss times, shown in 
Fig. The two time scales can indeed be comparable. 
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Figure 4. The distribution of time delays and arrival energies for the scenario "M82" 
for iron primaries, corresponding to Fig. [3 left panel. The horizontal line corresponds 
to the rectilinear propagation time. 
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Figure 5. The energy loss time as a function of energy for photo-disintegration on 
the combined cosmic microwave, infra-red and radio backgrounds. The solid line is 
for Helium nuclei, the dotted line for Carbon, the dashed line for Silicon, and the 
dash-dotted line for Iron. From Ref. |2M| . 

Furthermore, the spectra shown in Figs. El and El refer to the steady state situation 
relevant as long as the source emission characteristics do not change considerably on the 
time scale of typical time delays. This may cease to be realistic at low energies where 
time delays reach several hundred million years according to Fig. |3J Sources which 
were inactive at such times before the present would thus not be visible around 10 19 eV 
today, because such UHECR would not have reached us yet. This would produce a 
strong cut-off at low energies. 

We are here dealing with highly structured EGMF concentrated on scales of a 
few Mpc, with a typical r.m.s. strength B and coherence length L c . As long as the 
latter is not much smaller than the length scale over which the EGMF extends, the 
transition from rectilinear propagation to diffusion roughly occurs when the Larmor 
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radius r^{E) ~ E/(ZeB) becomes comparable to L c , or at 



For E < E t the time delay scales asymptotically as r(E) oc E 1 ^ 3 , whereas in 



the rectilinear regime t(E) oc E~ 2 [22J. According to Eq. (fl|). for light nuclei and 
L c ~ 100 kpc this transition occurs indeed between 10 19 eV and 10 20 eV, consistent with 
what is seen in Fig. H] for the branch with large delay times. This is consistent with 
the simulations of Ref. jTHj which has weaker and less extended fields than used in 
our simulations, but where in typical galaxy clusters protons diffuse nevertheless up to 
several 10 19 eV, as discussed in Ref. j3*T| . 

Another conspicuous feature in case of iron primaries is a rise of the spectrum 
towards lower energies at E ~ 4 — 6 x 10 19 eV, see Fig. EJ whose strength depends, 
however, considerably on the detailed EGMF realization. This energy roughly coincides 
with E inax /A where A = 56 for iron. It is thus given by the maximal energies of nucleons 
produced as secondaries by iron photo-disintegration. This is confirmed by the average 
atomic mass as a function of energy shown in Fig. 13 which shows a sharp drop at the 
same energy E ~ 6 x 10 19 eV due to these nucleons. Since such nucleons at energy E 
are produced by the flux of nuclei at energy AE, their effect is visible only for injection 
spectra not much softer than E~ 2 . The nucleon threshold is visible also in a branch of 
events with delay times r ~ 10 12 s around E ~ 60 EeV in Fig. 0] This branch is an 
almost identical copy of the analogous heavy nuclei branch with the same delay times 
at energies ~ 56 times higher. It is thus due to photo-disintegration of these nuclei 
into nucleons away from the source where the weak magnetic field hardly adds any time 
delay. 

Finally, again for iron primaries, the spectrum is considerably hardened in two 
energy intervals stretching over about half an order of magnitude, at E < 6 x 10 19 eV 
and at E < E max , see Fig. left panel. These two intervals are off-set from each 
other by roughly the charge of the iron primaries, Z = 26. A similar hardening is seen 
for proton primaries in the observed spectrum below ~ 3 x 10 20 eV, see Fig. 01 The 
hardening in these energy intervals can be interpreted by diffusive effects in terms of 
approximate analytical terms as follows: 

In the homogeneous case the diffusion - energy loss equation characterized by an 
energy dependent diffusion coefficient D(E) and (continuous) energy losses, dE/dt = 
-E/n OBB (E), has the solution [H2l ESI E2] 



for the flux j(E). Here, Q(E) is the injection spectrum of the discrete source at distance 
d and 




4\(E,E') 2 



(2) 
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Figure 6. Average atomic mass of observed UHECR as a function of energy predicted 
by scenario "M82" for iron primaries. The blue curves are for the EGMF surrounding 
the source as in Fig.^ left panel, and the red curve is without EGMF. The two panels 
correspond exactly to the two panels of Fig. [21 and demonstrate the dependence on 
source position. 



denotes an effective path length against energy losses. The diffusion coefficient scales 
as D(E) oc E a with a < 1 [22], and T\ oss (E) grows quickly with decreasing energy, see 
Fig. 03 As a result, X(E, E 1 ) increases with decreasing E. This gives rise to a change of 
slope in the flux Eq. (j2J), at an energy E d where X(E,E') ~ d: At energies E > E d 
where X(E, E') < d, the flux starts to become exponentially suppressed according 
to Eq. (j2J) which corresponds to the usual GZK-type cutoff. At energies E < E d 
the spectrum is modified by T\ oss (E) / X(E , E') 3 oc T\ oss (Ey l / 2 D(E)~ 3 / 2 relative to the 
injection spectrum. This factor is the ratio of the accumulation time T\ oss (E) and the 
volume X(E,E') 3 over which UHECR of energy E are distributed. When multiplied 
with the injection spectrum Q(E') in Eq. (J2J), it thus gives the UHECR density. In 
an EGMF with uniform statistical properties, D(E)~ 3 / 2 oc E~ 3a / 2 can increase faster 
with decreasing energy E than ri oss (.E)~ 1//2 decreases. This can lead to a steepening of 
the spectrum as observed, e.g., in the simulations of Refs. |221 ESI Oil an d discussed 
analytically in Ref. [33]. in contrast, in our inhomogeneous scenarios the effective 
diffusion coefficient D(E) depends more weakly on E because lower energy particles 
diffuse to larger distances from the source during their larger energy loss time. The 
correspondingly smaller EGMF partly compensates the decrease of D(E) with E. As 
a consequence, the factor T\ oss (E)~ 1 / 2 D(E) in general leads to suppression with 
decreasing energy. 

Since according to Fig. [3] the minimal T\ oss (E) is comparable to the source distance 
d in our scenarios, and since the transition from diffusion to straight line propagation 
occurs when D(E) ~ r^E) ~ L c which, as d, is of order Mpc, the energy E d is also 
roughly equal to the energy E t defined in Eq. This provides a rough explanation 
of the transition from diffusive hardening at low energies to softening at higher energies 
seen in the spectra of Figs. |2l Fig. [3J and Fig. [S] below. It should be kept in mind, 
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however, that our scenarios are highly inhomogeneous, the diffusion approximation 
strictly speaking does not apply far from the source, and these arguments can thus 
be considered only as approximate at best. Monte Carlo simulations are indispensable 
in this situation. 

In order to test how much the results depend on the detailed realization of EGMF 
and source, we also performed simulations in which the exact source position is varied 
within 1-2 Mpc within the general structures shown in Fig. ^ We find that the general 
tendency to spectral hardening at energies where UHECR diffuse in the vicinity of 
the source is always seen. However, the detailed shape of the predicted observable 
spectra depends considerably on the concrete structure of the EGMF. This is due to 
magnetic lensing [33] which leads to flux enhancements along specific propagation paths 
at specific energies, especially if there are field lines connecting source and observer. 
This is analogous to the direction dependence of the simulated nucleon spectrum from 
a source at the center of a magnetic slab obtained in Ref. [24J. Due to the unknown 
magnetic field structure these spectral details are in general unpredictable. This is 
demonstrated by comparing the two panels in Fig. |2 for iron primaries. The same is 
true for the detailed angular distributions. 
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Figure 7. The cumulative distribution of arrival direction off-sets from the source 
direction for UHECR above 4 x 10 19 eV in scenario "M82" for iron primaries, 
corresponding to Fig. [21 left panel (blue curve) and proton primaries, corresponding 
to Fig. |3| (red curve). 



Fig.|Ulshows the average atomic mass of observed UHECR as a function of energy for 
iron primaries. The nucleon spike towards low (A) at E ~ 6 x 10 19 eV is clearly visible, 
as discussed above. Below that energy there appears, however, also a relatively heavy 
component which survived excessive photo-disintegration. At all energies, the EGMF 
reduces the average atomic number due to increased photo-disintegration compared to 
undeflected propagation. Apart from this and the location of the nucleon spike, the 
detailed mass distribution is, however, quite sensitive to the EGMF structure, as seen 
by comparing the two panels in Fig. El 

Fig. El shows the cumulative distribution of arrival direction off-sets from the source 
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direction for the same "M82" scenario. Clearly, deflections are larger for iron primaries, 
whereas for proton primaries they are only a few degrees. 



Table 1. Approximate source emission power above ~ 10 19 eV in erg/s corresponding 
to the normalizations of the observed spectra in Figs.|3 left panel, |3 and|SJ 
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Figure 8. Same as Figs. [21 and [3] but for scenario "Virgo", shown in the right panel 
of Fig. ^ The left and right panels are for iron and proton primaries, respectively. All 
fluxes have been normalized at 60 EeV. 
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Figure 9. Same as Fig. but for scenario "Virgo" . 

We can also extract from the simulations the necessary source power corresponding 
to a particular normalization of the observed spectrum. For the normalizations shown 
in Figs. |2l El an d |S] below the source power is given in Tab. [T] The required power is 
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40. < E/EeV < 4000.: a= 31. 6± 23. 8± 23.8 
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Figure 10. Same as Fig. but for scenario "Virgo". This corresponds to the cases 
shown in Fig. [S] 

not strongly changed by the presence of the EGMF because we have normalized at the 
energy where the energy flux, which is proportional to the y-axis in Figs. El El and El is 
maximal. Note that for a given normalization of the observed spectrum, iron primaries 
require a somewhat higher injection power because a higher fraction of the injected 
energy is degraded below 10 19 eV during propagation due to spallation. For the case 
of "M82" most of the difference in required injection power between proton and iron 
primaries is due to the difference in normalization of the observed spectrum in these 
two cases, compare Figs. El and El Finally, the unknown EGMF structure translates 
into uncertainties in the required power of about an order of magnitude. Given these 
uncertainties of the injection spectrum and its continuation below 10 19 eV, these values 
are consistent with the UHECR emission power of active galaxies discussed in this 
context j22j, such as for M87, the main radio galaxy in the Virgo cluster |36j . 

Figs. El El and HH show spectra, atomic mass distribution, and the cumulative 
deflection angle distribution, respectively, for the "Virgo" scenario. For direct 
comparison with the "M82" scenario, the same injection parameters have been assumed. 
The effects of the EGMF surrounding the source, although more extended than in the 
scenario "M82", are here relatively smaller. This is because the delay time relative to 
the straight line propagation time now only reaches a maximum of ~ 10 at ~ 10 19 eV, 
compared to ~ 100 in the scenario "M82". The EGMF effects are, however, still 
significant and show the same generic features as discussed for the case "M82". In 
fact, for the same injection spectrum, at low energies, the suppression of the average 
atomic mass due to increased photo-disintegration by propagation in the EGMF is more 
severe than for the more nearby source "M82" , compare Figs. El and El This is because 
due to the larger propagation distance, fewer nuclei survive spallation. 

Note that the conventional pile-up in the spectrum around 100 EeV is clearly visible 
in the case of proton primaries without EGMF, see Figs. El and El right panels. 

Fig. ^2 shows an actual sky plot of arrival directions of events above 4 x 10 19 eV 
and above 10 20 eV for proton primaries in scenario "Virgo" . Several distinct images are 
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Figure 11. Sky distribution of UHECR arrival directions above 4 x 10 19 eV (left 
panel) and above 10 20 eV (right panel) for proton primaries in scenario "Virgo" . The 
source position is marked by an asterisk. The angular resolution was assumed to be 
1°. 



clearly seen which suggests significant magnetic lensing in the EGMF structure from 
Fig. Q right panel. This is not surprising given the fact that lensing in general sets in 
at energies higher than the transition energy to diffusion which, according to Eq. is 
roughly around 10 20 eV. The number of images and the deflection angles decrease with 
increasing energy, but are still significant at 10 20 eV and can be resolved with sufficient 
statistics. 

4. Conclusions 

We investigated the impact on cosmic ray observations above 10 19 eV of Mpc-scale 
magnetic fields of ~ 10~ 7 G strength surrounding ultra-high energy cosmic ray sources 
for the likely case that magnetic fields within a few Mpc around Earth are insignificant. 
We find that such source fields can strongly modify spectra and composition at Earth, 
especially for nearby sources for which the fields can considerably modify propagation 
times relative to both energy loss and photo-disintegration time scales and to the 
undeflected propagation time. We found the following generic features: 

The spectra are considerably hardened relative to the injection spectrum at energies 
below the usual GZK-like cutoff where energy loss distance and source distance become 
comparable. This is caused by an interplay between diffusion and energy loss: The 
flux of low energy particles is suppressed because diffusion spreads them out over a 
larger volume due to their much larger energy loss times. This is in contrast to the case 
of uniformly distributed magnetic fields which in general lead to a steepening of the 
cosmic ray flux below the GZK cutoff. A hardened sub-GZK spectrum from individual 
sources would be consistent with hints of a hard clustered component in the AGASA 
data between 10 19 eV and 10 20 eV jSj. 
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Furthermore, for a nucleus of atomic mass A as injected primary, due to the 
kinematics of the photo-disintegration reactions a nucleon peak appears at energy 
~ E max /A, where E mauX is the maximal nucleus injection energy. This effect is the 
more prominent the harder the injection spectrum. We also found that the details of 
spectra and composition depend significantly on the unknown details of the magnetic 
fields and the position of the source therein and can thus not be predicted. 

Further cutoffs towards low energies can be induced if the source is active only since 
a time smaller than the typical delay time at this energy. Obviously, the characteristics 
of such features also depend on unknown details of the source. 

Our simulations finally show that even for iron primaries, extra-galactic magnetic 
fields from large scale structure simulations are not strong and extended enough to 
explain the observed large scale isotropy of ultra high energy cosmic ray arrival directions 
in terms of a single nearby source. This would require more homogeneous fields such as 
in Ref. [2H1 - 

Next generation experiments such as the Pierre Auger Observatories [3] and 
the EUSO project jlj will accumulate sufficient statistics to establish spectra and 
distributions of composition and arrival directions from individual sources. A potentially 
strong influence of magnetic fields surrounding individual sources should thus be kept 
in mind when interpreting data from these experiments. This is true even if ultra-high 
energy cosmic rays arrive within a few degrees from the source position. 
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